Momentum transfer to small particles by aloof electron beams 
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The force exerted on nanoparticles and atomic clusters by fast passing electrons like those em- 
ployed in transmission electron microscopes are calculated and integrated over time to yield the 
momentum transferred from the electrons to the particles. Numerical results are offered for metallic 
and dielectric particles of different sizes (0-500 nm in diameter) as well as for carbon nanoclusters. 
Results for both linear and angular momentum transfers are presented. For the electron beam 
currents commonly employed in electron microscopes, the time-averaged forces are shown to be 
comparable in magnitude to laser-induced forces in optical tweezers. This opens up the possibility 
to study optically-trapped particles inside transmission electron microscopes. 

PACS numbers: 73.20.Mf,33.80.Ps,42.50.Vk,78.67.Bf 



I. INTRODUCTION 

Electromagnetic forces in optical tweezers are currently 
employed to trap small particles ranging in size from 
nanometers to several microns [l|, 3, and to manipu- 
late them in all spatial directions 3, This type of 
forces is also used to characterize the elastic properties 
of deformable tiny objects (e.g., living cells Q), to ob- 
tain quantitative information on mechanical properties 
at small length scales @, and in general, to fix the posi- 
tion of those particles so that they can be manipulated 
at will. 

In this context, transmission electron microscopy offers 
a potentially useful tool to study optically trapped par- 
ticles, providing excellent spatial resolution (sometimes 
below 1 A) when sub-nanometer electron beams are em- 
ployed [g], while allowing spectroscopic characterization 
with sub-eV accuracy. Actually, transmission electron 
microscopes are routinely exploited to probe local op- 
tical response properties Q, and more recently, also to 
determine photonic structures of complex materials Q . 

A major problem that may arise when combining elec- 
tron microscopy with optical tweezers or other types of 
optical trapping (e.g., optical lattices [U US, HH ) is that 
the passing electrons can kick the particles out of the 
trapping locations (see Fig. [I}. In this work, we show 
that the momentum transferred from the passing elec- 
trons to the particles can be well below the threshold 
needed to kick them out for commonly employed trapping 
laser intensities, although a detailed comparison between 
trapping forces and electron-induced forces suggests that 
both weak and strong perturbation regimes are possible 
depending on the distance between the particles and the 
beam, all of them within the range that allows a suf- 
ficiently large electron-particle interaction as to perform 
electron energy loss spectroscopy (EELS) with significant 
statistics for in vacuo optically-trapped particles. 

The moving electrons can be in fact regarded as a 
source of evanescent electromagnetic field that probes the 
sample locally, and in this sense, they can be also used 
to produce deformation in elastic particles, oscillations 
of trapped particles around their equilibrium positions, 



FIG. 1: (color online). Schematic representation of the pro- 
cess considered in this work: a fast electron moving with im- 
pact parameter b and velocity v with respect to a polarizable 
particle transfers momentum Ap = (Ap x , Ap z ) to the particle 
via electromagnetic interaction. 



and other interesting effects associated to the transfer of 
momentum within accurately controlled spatial regions. 

The present work addresses, in a quantitative way, the 
issue of both linear and angular momentum transfer from 
an electron beam to small particles. This applies to opti- 
cally trapped particles, as mentioned above, but also to 
other forms of trapping, like particles deposited on solid 
substrates, or particles trapped by a tip (e.g., in a STM 
set up). 



II. THEORY 

The electromagnetic force exerted on a particle in vac- 
uum is given by the integral of Maxwell's stress tensor 
over a surface S embedding the particle [l2j as 



F(i) 



— / ds[E(s,t)E(s,t) • n + H(s, t) H(s, t) • n 
47r Js 



|E(s,t)| 2 + |H(s,0| 2 )], 
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where n is the surface normal and Gaussian units are 
used. The momentum transferred to the particle, Ap, is 
obtained by integrating of F(t) over the time. This yields 



Ap = J F(t) dt = F M 



did, 



(1) 



where 



F(w) = T^Re{/ dBtE( B| w)(E(s,w)-A)* (2) 
+H(s,w) (H(s,w)* • n)* 
-|(|E(s, W )| 2 + |H(s, W )| 2 )]}, 

and the Fourier transform is defined as E(r, u>) — 
J dtE(r,t)exp{iut}. 

The force acting on the particle is due in part to radi- 
ation emitted as a result of interaction with the electron 
and in part to the reaction force experienced by the pro- 
jectile. For small particles, the effect of radiation emis- 
sion is negligible and the trajectory is deflected by an 
angle w Ap/mv, where m and v are the mass and veloc- 
ity of the electron. Non-retarded calculations have shown 
that this angle is too small to be easily measured [T3 |. 



A. Small particles 

Let us first consider a small isotropic particle suffi- 
ciently far away from the electron beam as to neglect 
higher multipoles beyond induced dipoles. The particle 
is then characterized by its frequency-dependent polar- 
izability a(u>), and the force exerted by each frequency 
component of the external field E(r, ui) reduces to [lij 

F(w) = Re{a Y, E T *(r,w)V[£f ( 3 ) 



This expression can be derived from Eq. |(5J) by consider- 
ing an integration surface arbitrarily close to the object 
and by using the expressions for the electric and magnetic 
fields induced by a small polarizable particle in terms of 
its polarizability a. For an electron moving with velocity 
v towards to positive z direction and passing by the ori- 
gin at t = 0, the external field is readily calculated from 
Maxwell's equations to yield 



E ext (r,w) 



-2euj : 



V 2 "f 



luR^R i 

vj R 7 



K (— )i],(4) 
«7 



where R = (x, y) and 7 = 1/yT — v 2 /c 2 . Inserting Eq. 
dH) into Eq. f3}, one obtains 

F M = Sh^M * + 27lm{a} /(^) z],(5) 



where 



/(C)=*r 2 (C) + ^ 2 (C)/7 2 



and the particle is taken to be situated at R = (—b, 0) 
with respect to the beam (see Fig.[TJ). 

Symmetry considerations lead to the conclusion that 
Rayleigh scattering of the external-electron evanescent 
field ([4]) produces a radiation pattern with inversion sym- 
metry with respect to a plane perpendicular to the tra- 
jectory. This means that the overall transfer of momen- 
tum to the induced radiation is zero in the small-particle 
limit, so that Ap z accounts for all momentum transfer 
to the moving electron along z. Then, the contribution 
of each to component to the electron energy loss rate is, 
within the non-recoil approximation valid for sufficiently 
energetic electrons, vF z (u). Actually, one finds that the 
identity vF z (u) = KujP(uj) is satisfied, where P(u>) is 
the frequency-resolved loss probability as previously ob- 
tained for small particles [151 ]. As a consequence, F z van- 
ishes in the u> — > limit, since P(w) remains finite. 

This behavior is quite different from F x , which goes 
to a finite value for small w's, namely F x (u> = 0) = 
4e 2 Re{a(0)}/i> 2 & 3 . (Incidentally, momentum transfer 
along x produces negligible energy transfer in the non- 
recoil approximation.) This latter formula can be used 
to derive a close expression for Ap x valid for arbitrarily- 
large, finite objects in the large impact parameter limit. 
In that case, only small w's contribute to F(w), due to 
the effective exponential cut-off imposed by the modi- 
fied Bessel functions Kq and K\. This means that only 
long wavelengths are relevant (to which the object ap- 
pears as small), so that it can be described by its static 
polarizability. Then, the u> integral can be performed 
numerically to yield 

frr „„ 1.85055, e 2 Re{a(0)} 
Ap x = (5.55165 7 + ) TT^ 1 - (6) 

7 vr 

For comparison, the momentum transferred to a charge 
e at a distance b from the beam is Ap = —(2e 2 /bv)'k. 

The large-6 limit given by Eq. ((6]) is compared in Fig. 
[2] with more detailed calculations that include higher- 
multipole moments, as described below. Also, the small 
particle limit of Eq. ([5]) is discussed in Fig. [31 



B. Arbitrary size 

For larger particles or for close electron-particle en- 
counters, higher multipoles become relevant in the in- 
duced forces [161 ]. Then, it is convenient to express the 
evanescent field of the electron in terms of multipoles 
centered at the particle, so that the external electric 
and magnetic fields admit the following decomposition 

MM- 

E cxt (r,c) = ]>>f <° xt L - Ul ' cxt V x L]j L (kr) 

L K 

and 

H cxt (r,c) = -]>>f cxt L + A^f '° Xt V x L]j L (kr), 
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FIG. 2: (color online). Momentum transfer to small spherical 
particles by a passing 200-keV electron as a function of the 
distance from the trajectory to the center of the spheres b. 
The momentum transfer has been scaled using the velocity 
v — 0.7c, the sphere radius a, and the impact parameter 
b. The perpendicular component of the momentum transfer 
with respect to the trajectory Ap x (solid curves) has been 
represented for spheres of radius a — 10 nm, 50 nm, 200 nm, 
and 500 nm (notice the rapid increase in Ap x near b = a). The 
parallel component Ap z (dashed curves) is only shown for a — 
200 nm and 500 nm. Dielectric alumina spheres and metallic 
silver spheres are considered (left and right plot, respectively), 
respectively. The large b limit for perpendicular momentum 
transfer [Eq. ©] is shown by horizontal dotted lines. 
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and 



C± = y/(l±m + l)(lTm). 

1 if to > 0, s m = (-!) m if m < 0, and C$ is 
the Gegenbauer polynomial [18j |. The impact parameter 
b is defined in Fig. [TJ 

The induced field around the particle is given by simi- 
lar expressions obtained by substituting ijj v £ ext by new co- 
efficients ^ lnd , and ji by the Hankel function 
In particular, L = (I, to) is conserved for spherical par- 
ticles and one has a linear dependence V , L ind = titp^ ext , 
where t" are scattering matrices that are given by ana- 
lytical expressions in the case of ho mog eneous particles 
of dielectric function e and radius a [151 ] : 
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FIG. 3: (color online). Particle size dependence of the mo- 
mentum transfer normalized to the particle mass M under 
the same conditions as in Fig. [2] small particle limit (dashed 
curves) versus full multipole calculation (solid curves). The 
particle is made of AI2O3 (density p — 4.02 g/cm 3 ), the elec- 
tron energy is 200 keV, and the distance from the trajectory 
to the particle surface is 10 nm. Dotted curves show the mo- 
mentum transferred from light in an optical trap (see text for 
details) . 



where L — (l,m), k — ui/c, j'i(fcr) = rj';(fcr)Yi,(r), L = 
— ihr x V is the orbital angular momentum operator, and 
?/>r' ext (for v — E, M) are multipole coefficients given by 

Gam 



, M,cxt 
i E.cxt 



1(1 + 1) He 



S1/7 



,M 



-MPo)pd'i(Pi) + Poj'i(po)ji(pi 



,(+> 



(Pa)Pij'i(Pi) - Po[h\ + \po)]'ji(p 



and 



t E _ -ji{po)\piji{pi)]' + 4po.n(po)Yji(pi) 

1 hf +) 0*)[piii(pi)]' - e[p h ( l +) (po)}'ji(piy 

where po — ka, p\ — po\/e. with Imjpi} > 0, and the 
prime denotes differentiation with respect to po and p\. 

At this point, it is convenient to write the operators L 
and (l/fc)V in matrix form. One finds 

L?L = ^2 L LL'jL' 

V 



and 



respectively, where 

L L L' = hSi t i'[C + ^ m+ i, m /(x- iy)/2 

+C_ (5 m _i, m '(x 4 iy)/2 + mi5 m , m / zj, 



z-N 



LL' 



■A IK j- A , [V + m){V -m) 
i<W(W + *-i,J') {2l , _ 1)(2Z , + x) , (8) 



and the x and y components of N are obtained from ([5]) 
by rotating the reference frame using rotation matrices 



4 



for spherical harmonics [19j . Exactly the same matri- 
ces as above apply to L and (l/fc)V acting on Hankcl 

functions . Furthermore, these matrices satisfy the 
properties L + = L and N + = — N. 

Now, the electric field admits an expansion of the form 

E cxt (r, W )=^ErjL(fcr), 

L 

where the coefficients 

E r = y, LLL^< oxt + i £ N^ L „ x \,i„ L ,^\ 

V L'L" 

are obtained from the above expressions. Similar formu- 
las are obtained for H cxt and for the induced fields E md 
and H md in terms of multipole coefficients. Finally, we 
insert them into Eq. ((2]) and perform the integral over a 
sphere in the s — > oo limit. Then, the first two terms in- 
side the integrand give a vanishing contribution because 
the induced far-field is transverse. The remaining part of 
the integral can be recast, noticing that only real terms 
must be retained, 



F(w) 
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JJCXt 
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(Hi nd n(i-(-i) ( ) 

(HS?)*](l-(-l) ,# ) 
(H£ d )*])}, 



where 



h LL , = J dnY*,{Q)h(n)Y L {Q) 



(10) 



andn(fi) = [(x+iy) Yx-x/y/2- (x-iy) Y u /y/2+ 

zYio] is the radial vector as a function of the polar di- 
rection f2. 



III. RESULTS AND DISCUSSION 

Fig. [5] shows the dependence of the momentum trans- 
fer on electron impact parameter b for alumina and sil- 
ver spheres of different radius, as calculated from Eqs. (Q]) 
and © ■ Measured optical data have been used for the di- 
electric function of these materials [?J • One can observe a 
nearly exponential decay of the momentum transfer with 
b. Besides, the momentum transferred along the direction 
of the electron velocity vector (Ap z , dashed curves) is 
generally smaller than the remaining perpendicular com- 
ponent (Ap x , solid curves), which finds an explanation 
in the fact that the contribution of these components to 
the energy loss Tiuj is vAp z + (Ap x ) 2 /to, where m is the 
electron mass: since mv 3> Ap, Ap x is allowed to take 
larger values than Ap z for each fixed u). 

Notice also that Ap x converges quickly to the large b 
limit [Eq. [51 dotted curves] , producing a finite result un- 
der the scaling of Fig. [5J unlike Ap z , which goes faster to 
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FIG. 4: (color online). Momentum transferred from a 5-keV 
electron to a Ceo cluster as a function of impact parameter b. 
The momentum is normalized to the cluster mass M. 



for large b. In this limit, the electron induces a dipole in 
the particle directed towards the electron, which results 
in an attractive force between these two similar to the 
image potential at surfaces [20| , leading to a momentum 
transfer Ap w Ap x x. For small metallic particles and 
closer encounters this picture is no longer valid and Ap x 
can actually reverse its sign and have a net repulsive be- 
haviour (e.g., in Fig. [2] for Ag particles of radius a = 10 
nm and also for the fullerenes of Fig. SJ) . 

A more detailed analysis of the magnitude of the mo- 
mentum transfer effect is given in Fig. [3l The momentum 
transfer is normalized to the particle mass M and the re- 
sult is the change in the particle velocity induced by the 
passage of the electron as a function of particle radius a. 
The trajectory of the 200-keV electron under considera- 
tion passes 10 nm away from the surface of the spherical 
alumina particles. The full-multipole calculation [Eqs. 
([1]) and ([9]), solid curves] agrees well with the small par- 
ticle limit [Eqs. ([1]) and (|5|), dashed curves] when a is 
much smaller than b — a = 10 nm. Even though the 
electron-particle interaction increases with the radius a, 
the actual change in the particle velocity shows a nearly 
exponential decay with increasing a. 

In a situation where the particle is trapped by lasers 
(e.g., in optical tweezers [|[ or in optical stretchers [|[), 
one should compare the interaction with the electrons to 
the interaction with the laser light. To this end, we will 
consider a trapping cw-Ti:sapphire 100-mW laser emit- 
ting at a wavelength A = 785 nm and focused on a region 
of radius Rf = 10 /im. Furthermore, we will contem- 
plate the momentum transferred by the laser during the 
average time span At between two consecutive passing 
electrons in a transmission electron microscope operating 
at a current of 1 nA. The particle polarizability a is all 
that is needed to calculate light forces for the small radii 
under discussion (a <C A), according to Eq. (J3j) . Now, 
for real a this equation defines a conservative gradient 
force that responds to the potential — (a/2)|E| 2 , where 
E is the laser light field, whereas the imaginary part of a 
represents photon absor ptio n by the particle that trans- 
lates into light pressure [21| . These two components are 
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represented separately in Fig. [3J after multiplication by 
At/M (dotted curves). The light pressure contribution 
is calculated for an incidence plane wave with the same 
photon flux as the laser at its focus. The gradient force 
component is obtained from the maximum force in the fo- 
cus region assuming a Gaussian profile for the laser field 
intensity (i.e., |E| 2 oc exp[-R 2 /(R f / In 2) 2 ]). Finally it 
is convenient to define the polarizability from its relation 
to the scattering matrix, which upon inspection permits 
writing a = 3if /2fc 3 . Unlike the well-known expression 
[l2| a = a 3 (e — l)/(e + 2), the former relation predicts 
a non-zero value for Im{a} even for particles with real 
e (like our alumina spheres), arising as a pure retarda- 
tion correction associated to radiation scattering (this is 
actually the origin of the light pressure component of 
Fig. [3]) • (Incidentally, gravity would produce a velocity 
change gAt — 1.56 nm/s, which is well compensated for 
in currently available optical trapping systems.) 

An important conclusion that can be extracted from 
Fig.[3jis that the crossover of trapping light into the main 
source of momentum occurs for particles of 20 nm in di- 
ameter when the electrons pass at a distance of 10 nm 
from the particles surface, thus allowing one to perform 
energy loss analysis of the transmitted electrons with sig- 
nificant statistics. Therefore, transmission electron mi- 
croscopy can be combined with in- vacuo optical trapping 
to study particles of sizes above some tens nm. 

While the transfer of momentum by the trapping light 
occurs in a continuous smooth fashion, the electrons de- 
posit all of the momentum during a small time interval 
~ a/v (<C At = 0.16 ns for 1 nA electron current). How- 
ever, the change in particle velocity per electron (vertical 
scale in Fig. [3J produces a minute particle displacement 
during At (smaller than 1.6 x 10~ 9 nm -C a), and there- 
fore, the effect of the passing electrons is experienced by 
the particle as a nearly continuous source of momentum 
that is describable by an average force Ap/At. Actually, 
Fig. [3J suggests that using more intense electron beams 
(with even smaller impact parameters) acting during pe- 
riods of the order of one second will still not produce 
ejection of the particles from their trapping locations. 

It should be stressed that the momentum transfers that 
we have calculated using classical electromagnetic the- 
ory must be understood as the average value over many 
incoming electrons, since the actual strength of the in- 
teraction is not large enough as to guarantee that many 
photons are exchanged between each electron and a given 
particle. Like in aloof EELS experiments [g], most elec- 
trons will not interact with the particles at all, so that the 
present results must be understood under the perfectly 
valid perspective of an statistical average performed over 
many beam electrons. The quadratic deviation from 
these average forces can play also a role (similar to strag- 
gling in stopping power theory), but this subject is left 
for future consideration. 

We have also studied momentum transfer to Cgo clus- 
ters (Fig. 2]). The scattering matrices t\ have been ob- 
tained within the discrete-dipole approximation (22l . |23| , 
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FIG. 5: (color online). Change in the particle angular velocity 
as a result of the torque exerted by the electron. The particle 
is made of AI2O3 (density p = 4.02 g/cm 3 ), the electron en- 
ergy is 200 keV, and the distance from the trajectory to the 
particle surface is 10 nm. 



where each carbon atom is described by an induced dipolc 
whose polarizability is fitted to reproduce correctly the 
measured optical response of graphite Further de- 
tails concerning the procedure followed to obtain t\ will 
be given elsewhere [24[. At relatively small interaction 
distances 6, the z component of the momentum is larger 
than the x component and the latter is negative. These 



are effects that can be hardly found in the above exam- 
ples and that originate in high-order multipoles (actually, 
I < 5 are needed for convergence within the range of b 
under consideration). Even at a distance of 9 nm (no- 
tice that Ceo has a diameter of only 0.7 nm) the change 
in velocity produced by the passing electron can be sub- 
stantial. Therefore, the interaction of fast electrons with 
small clusters can produce dramatic effects if these are 
not mightily bound by a mechanism stronger than optical 
trapping. 

Finally, the passing electron can induce a torque on 
the particle that changes its angular momentum (AL y ) 
and makes it rotate. This is the effect discussed in Fig. 
[5] Like the electromagnetic force above, the torque G 
is obtained from the integral of Maxwell's stress tensor 
[l2l ] , which yields 

G = Re{ / dsUh ■ E)*(E x n) + (n • H)*(H x n)]}, 

4?r Js 

where s is the radius of a sphere where the particle is 
embedded. Proceeding in a way similar to the derivation 
of the force presented in Sec. |TIJ this expression can be 
written in terms of multipole components as 



(11) 



Imm' 



; \/(l + m + l)(l — m)6 m+li7n i (x - iy) 



+ -^(l-m + i)(l + rn)5 m - ltm >(x + iy) 
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x [^ ind (C' nd )*+^ ind (^ nd )* 
+i^f d (^ Xt )* + i^ ind (^ xt )1}. 

Eq. (fTTjl has been used to obtain Fig. [5j which shows 
the change in angular velocity of the particle per inci- 
dent electron, Aft = AL y /I, where / = (2/3)a 2 M is the 
moment of inertia of the alumina sphere. 

Averaging over the electrons of a 1 nA electron beam 
passing at 10 nm from the surface of an alumina sphere 
of radius a = 20 nm, one finds an angular acceleration 
of 39 MHz/s. Under these conditions, the linear mo- 
mentum transferred by the electrons can be absorbed by 
the trapping light, as discussed above. However, the an- 
gular momentum is not absorbed, and the particle will 
spin with increasing angular velocity until either the cen- 
trifugal force breaks it apart or radiation emission at the 
rotation frequency (vacuum friction) compensates for the 
electron-induced torque. 

In conclusion, we have shown that fast electrons follow- 
ing aloof trajectories (i.e., without direct overlap with 
the sample) in a transmission electron microscope can 
exert time-averaged forces on small particles of similar 
magnitude as those forces associated to trapping in opti- 
cal tweezers and stretchers, and therefore, this effect can 
be used for analytical studies of mechanical properties 



of such particles, while electron energy loss spectra can 
be actually taken without causing ejection of the par- 
ticles from their trapping positions. The experimental 
challenge consists in combining electron microscopy with 
optical trapping in vacuum, and we hope that this work 
can contribute to stimulate further research in this direc- 
tion. Besides, the linear and angular momentum trans- 
fers discussed here can play a relevant role in other types 
of trapping, like for particles deposited on a substrate, 
or for particles trapped by a tip, where the forces and 
torques induced by the passing electrons can push, pull, 
reposition, or reorient the trapped particles. 
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